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Abstract 

The near threshold expansion of generalized sunset-type (water melon) diagrams 
with arbitrary masses is constructed by using a configuration space technique. We 
present analytical expressions for the expansion of the spectral density near thresh- 
old and compare it with the exact expression obtained earlier using the method of 
the Hankel transform. We formulate a generalized threshold expansion with partial 
resummation of the small mass corrections for the strongly asymmetric case where 
one particle in the intermediate state is much lighter than the others. 



1 Introduction 



Perturbation theory remains the most accurate and reliable tool of theoretical analy- 
sis in elementary particle phenomenology. Even the most promising approach based on 
the direct calculation of physical observables on the lattice still cannot compete with 
perturbation theory supplemented with some non-perturbative additions as condensates 
or renormalons. The accuracy of experimental data for tests of the Standard Model 
permanently improves and demands the improvement of the accuracy of theoretical pre- 
dictions 0. Within perturbation theory this demand implies the calculation of Feynman 
diagrams with an increasing number of loops (see e.g. 0). At present, Feynman diagrams 
of rather general form can be obtained analytically at most up to three loops. Any of 
these calculations is still unique and involves sophisticated algorithms and an extensive 
usage of computer resources (see e.g. ||, U). Even the notion of the analytical calcula- 
tion is changing now because the analytical results are expressed through rather special 
transcendental numbers which eventually have to be calculated numerically. The striking 
example is the calculation of massive three-loop bubbles ||] where the analytical results 
for the diagrams are represented through the quantity which is given by a two-fold infinite 
sum (see also f|). Not all of the multiloop Feynman diagrams are equally complicated. 
There are some topologies of n-loop diagrams or special kinematic regimes for the ex- 
ternal momenta which are much simpler to compute than the general n-loop case (see 
e.g. 0). The increasing complication of loop calculations within perturbation theory is 
more and more circumvented by using different kinds of expansion for the cases where the 
kinematics allows one to find a small parameter. Heavy Quark Effective Theory (HQET) 
and threshold expansions for heavy particle production are recent examples for such a 
situation. These expansions are mostly used for phenomenological applications where it 
even suffices to give approximate numerical values for the unexpanded multiloop integrals. 
In some cases one can obtain an analytical expansion that can be used for the simplifi- 
cation of the calculations or for the evaluation of the exact parameters of the effective 
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Lagrangians through matching conditions. 

In this paper we discuss a method to calculate the near threshold expansion of the 
spectral density for the class of multiloop sunset-type Feynman diagrams termed water 
melon diagrams J7|. These diagrams has recently drawn some attention as a laboratory 
for testing some advanced methods in loop calculations || However, water melon 
diagrams also have numerous phenomenological applications. The most important ap- 
plication of these diagrams (as well as spectacle diagrams as their first generalization) 
is the calculation of the effective potential in quantum field theory both at finite tem- 
perature and at zero temperature [10, 11, 12, [13], 14, 15]. Among other applications one 
can name the sum rule analysis of baryons in QCD both in the massless approximation 
and with finite mass (the leading order correlator being the standard sunset diagram) 
D], 18|, [19| or the treatment of baryons in the large iV c limit of QCD |2^, ^TJ which 



requires the calculation of (N c — l)-loop water melon diagrams in the leading order. An 
interesting application is the calculation of dibaryon properties in operator product ex- 
pansion within the sum rule approach which is important for understanding many-quark 
bags in the nucleus E2L E3. The water melon diagrams emerge in chiral perturbation 



theory |24], pqj . The corresponding integrals appear for the correlator of the effective 
operators related to the mixing of neutral mesons in flavour dynamics |26[] and in the 
sum rule analysis of hybrid mesons |27| . Water melon diagrams constitute an important 
part of the contributions generated by the recurrence relations of the integration-by-part 
techniques for three-loop diagrams |28j (for the recent progress see e.g. Ref. p9|). It is a 
rather hot topic and a dynamically developing area of loop calculations where new results 
are frequently reported (see e.g. Refs. |3(| |3l], [32f ). 

In our opinion, the method presented in Ref. |7j completely solves the problem of 
computing this class of diagrams. The method is simple and reduces the calculation of 
a n-loop water melon diagram to a one- dimensional integral which includes only well- 
known special functions in the integrand (Bessel functions of different kinds) for any 
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values of the internal masses. The technique is universal. The most interesting part of 
our analysis of water melon diagrams is the construction of the spectral decomposition 
of water melon diagrams, i.e. the determination of the discontinuity across the physical 
cut in the complex plane of the squared momentum. In this context a novel technique 
for the direct construction of the spectral density of water melon diagrams based on an 
integral transform in configuration space was presented in Ref. 0. In the present paper 
we develop some explicit expansions for water melon diagrams and compare them with 
the exact results. The purpose of this paper is a three-fold one, namely: 

• to demonstrate the ease with which threshold expansions to arbitrary space-time 
dimensions and a general number of internal lines with arbitrary masses can be 
generated within a configuration space technique. The entire construction reduces 
to algebraic operations while the only integral encountered is a simple integral of 
the type of Euler's Gamma function. 

• to generate explicit forms of threshold expansions and to analyze their convergence 
properties. The explicit forms of threshold expansions for the simple sunset can be 
compared with threshold expansions that are obtained in momentum space. 

• to consider a case when one mass is much smaller than the others. This is an 
important generalization of the threshold expansion and can be analytically done 
within the configuration space technique. 

The present analysis of the threshold expansion of water melon diagrams shows the effec- 
tiveness of the configuration space technique for this topology class of Feynman diagrams. 
On the other hand the obtained explicit results can be useful for a variety of applications 
which include the evaluation of water melon type diagrams. The technique can readily 
be generalized to some other simple topologies of n-loop diagrams. 

The paper is organized as follows. In Sec. 2 we provide the tools necessary for the 
following calculations. Sec. 3 deals with the threshold expansion. Here we outline our 
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main strategy for this expansion in configuration space and in the subsections that follow 
we give examples for the sunset diagram and the water melon diagram with four or more 
internal lines. In Sec. 4 we concentrate on the strongly asymmetric case where one of 
the masses is much smaller than the others. We introduce the procedure of resummation 
of the contributions of this smallest mass and show explicitly how we get to a closed 
expression for the spectral density even in this case. We give examples for a strongly 
asymmetric mass arrangement for the water melon diagrams with two (degenerate case), 
three, and four or more internal lines, and compare the resummed results with the pure 
threshold expansion as well as with the exact spectral density. In Sec. 5 we discuss how 
to recover the non-analytic part of the polarization function through the spectral density 
near threshold. In Sec. 6 we give our conclusions. 



2 Basics about the configuration space technique 

We start with a brief outline of the technique which we use in this paper 0. The po- 
larization function U(x) of a water melon diagram including n internal lines with masses 
rrii, i G {1, . . . , n} in configuration space is given by the product 

n 

U(x) = l[D(x,m l ). (1) 

i=i 

The propagator D{x,m) of a massive line with mass m in .D- dimensional (Euclidean) 
space-time is given by 

1 r e ip ^d D p _ (mx) x K\(mx) 
(X ' m)= (2^pi p 2 + m 2 = (2tt) x + 1 x^ [ } 

where we write D = 2A + 2. K\(z) is a McDonald function (a modified Bessel function of 
the third kind, see e.g. In the limit m — > the propagator in Eq. (Q) simplifies to 

D(x0) - 1 [ e^d° P _ T(A) 
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where T(A) is Euler's Gamma function. Eq. ([I]) contains all information about the water 
melon diagrams and in this sense is the final result for the class of diagrams under con- 
sideration. It is of course known and was used since long ago f34}] . Of some particular 
interest is the spectral decomposition of the polarization function IT(a;) which is connected 
to the particle content of a given model. The spectrum of particles is contained in the 
function p(s) related to the polarization function through the dispersion representation 



The dispersion representation of the polarization function in configuration space reveals 
the analytic structure of the polarization function U(x). For applications, however, one 
may need the Fourier transform of the polarization function U(x) given by 



(We use the same notation for the function and its Fourier transform because we think 
that this will cause no confusion). Because of the Lorentz invariance of the propagator 
the angular integration in Eq. (|5|) can be done explicitly in D-dimensional space-time (for 
a generalization to tensor propagators see Ref. f7|). The result reads 



where J\(z) is the usual Bessel function and d D Q is the rotationally invariant measure on 
the unit sphere in the D-dimensional (Euclidean) space-time. Note that the polarization 
function U(x) as well as its Fourier transform 11(g) are only functions of the absolute 
value \x\ and |g|. The same is of course valid also for the other occurring functions. To 
simplify the notation we often write x — \x\ and q = \q\ for these absolute values. 

Note that propagators of particles with non-zero spin in configuration space repre- 
sentation can be obtained from the scalar propagator by differentiation with respect to 
the space-time point x. This does not change the functional x-structure and causes only 
minor modification of the basic technique (for details see Ref. and references therein). 




(4) 




(5) 
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Our final representation of the Fourier transform of a water melon diagram is given by 
the one-dimensional integral 

U(q) = 2tt a+1 J™ (^j X J x (qx)U(x)x 2X+1 dx (7) 

which is a special kind of integral transformation with a Bessel function as a kernel. This 



integral transformation is known as the Hankel transform [|35], 3(|. The representation 
given by Eq. (f7|) is quite universal regardless of whether tensor structures are added or 
particles with vanishing momenta are radiated from any of the internal lines. Note in this 
context that Bessel functions are objects well-studied during the last two centuries. They 
therefore can be added to the list of elementary functions. 

Because the dispersion representation of the polarization function in configuration 
space (or the spectral density of the corresponding polarization operator) has the form 
given in Eq. (|j), the analytic structure of the polarization function U(x) can be determined 
directly in configuration space without having to compute its Fourier transform first. The 
transformation in Eq. (|J) turns out to be a particular example of the Hankel transform, 



namely the i^-transform |35], [36fl . As pointed out in Ref. 0, the inverse f^-transform in 
this case is given by 

p(s) = / /A(Cv^)n(C)C A+1 rfC (8) 

IS 1 Jc—ioo 

where I\{z) is a modified Bessel function of the first kind and the integration runs along 
a vertical contour in the complex plane to the right of the right-most singularity of n(£). 
The inverse transform given by Eq. (|8|) completely solves the problem of determining the 
spectral density of the general class of water melon diagrams with any number of internal 
lines and different masses by reducing it to the computation of a one-dimensional integral. 
For n internal lines with equal masses m the spectral density reads explicitly 

rr>^ n rc+ioo 

*»> = i(2l , r -^n jL h(CV ~ s) (^("0>"C 1 - ( "- 1 ^C. (9) 

In contrast to our technique, in the standard, or momentum, representation the polariza- 
tion function n(g) is calculated from a (n — l)-loop diagram with (n — 1) D-dimensional 



integrations over the entangled loop momenta which makes the computation difficult when 
the number of internal lines becomes large. 



3 Threshold expansion 

With our method described in detail in Ref. J7j the s-dependence of the spectral density 
can be calculated by a one- fold numerical integration according to Eq. (||). The numerical 
integration in Eq. ([§D can be done for arbitrary space-time dimensions and a general 
number of lines with arbitrary masses. In this sense this is the most efficient representation 
for the spectral density of the water melon diagram. However, we can also develop an 
explicit expansion near the threshold with any desired accuracy. It does not require 
any complicated integrations at all. The corresponding expansion of Eq. (^|) can then 
be compared with series expansions near the production threshold obtained with the 
traditional momentum space technique. We stress that we are only interested in the 



spectral density because it is the main object for physical applications (see e.g. Refs. [j37| ). 
For practical reasons we start with Eq. (|7|). The polarization function 11(g) as written 
in Eq. (|[|) is, in general, UV divergent. The divergence can be subtracted by using the 
power series expansion of the weight function (qx/2)~ x J\(qx) to an appropriate order 
which will be added and subtracted to this weight. This leads to a g 2 -dependent power 
series of divergent subtraction terms plus an UV finite subtracted integral (see e.g. [|7|, j38f). 
But because the subtraction terms do not contribute to the spectral density, we can 
avoid this subtraction at all. In order that the formally written expressions make sense 
they are supposed to be dimensionally regularized. We use the simplified or unorthodox 
dimensional regularization method for water melon diagrams (see Ref. 0). 

The threshold region of a water melon diagram is determined by the condition q 2 + 
M 2 ~ where q is the Euclidean momentum and M = Y^,i m i is the threshold value for the 
spectral density. We introduce the Minkowskian momentum p defined by p 2 = —q 2 which 
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is an analytic continuation to the physical cut. Operationally this analytic continuation 
can be performed by replacing q — > ip. To analyze the region near the threshold we use the 
parameter A = M — p which takes complex values. The parameter A is more convenient 
in Euclidean domain while the parameter E = — A = p — Mis the actual energy counted 
from threshold which is used in phenomenological applications. The spectral density as a 
function of E is written as p{E) = p((M + E) 2 ) in the following. The analytic continuation 
of the Fourier transform in Eq. (|7|) to the Minkowskian domain has the form 

n(p) = 2tt a+1 J™ (^j X J x (ipx)Yi{x)x 2X+l dx. (10) 

For the threshold expansion we have to analyze the large x behaviour of the integrand. 
It is this region that saturates the integral in the limit p — > M or, equivalently, E — > 0. 
It is convenient to perform the analysis in a basis where the integrand has a simple large 
x behaviour. The most important part of the integrand is the Bessel function J\(ipx) 
which, however, contains both rising and falling branches at large x. It resembles the 
situation with elementary trigonometric functions sin (2) and cos(z) to which the Bessel 
function J\(z) is rather close (in a certain sense). Indeed, cos(z) (or sin(z)) is a linear 
combination of exponentials, namely 

cos(z) = \ (e" + e-«) (11) 

and has also both rising and falling branches at large pure imaginary argument: the ex- 
ponentials show simple asymptotic behaviour e ±z at z = ±zoo. The analogous statement 
is true for J\(z) which can be written as a sum of two Hankel functions, 

Jx(z)^ 1 -(H+(z) + H^z)) (12) 

where H x (z) = J\{z) ± iY\{z). The Hankel functions H x (z) show simple asymptotic 
behaviour at infinity, 

Hf{iz) ~ z~ 1/2 e ±z . (13) 
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Accordingly we split up II(p) into II(p) = IT + (p) + II (p) with 

n±(p) = v X+1 J™ (jy) X H^{tpx)U{x)x 2X+1 dx. (14) 

The two parts IT^p) of the polarization function IT(p) have completely different behaviour 
near threshold which allows one to analyze them independently This observation makes 
the subsequent analysis straightforward. We first consider the contribution of the part 
Il + (p). The behaviour at large x is given by the asymptotic form of the functions which 
we simply write up to the leading terms as 

H + (ipx) = </— -e-^fl + Oix" 1 )), K(mx) = J— e - mx (l + O^" 1 )). (15) 
V mpx V 2mx 

The large x range of the integral (above a reasonably large cutoff parameter A) has the 
general form 

n+(M - A) ~ / x- a e- i2M - A)x dx (16) 

J A 

where 

o=(n-l)(A + l/2). (17) 

The right hand side of Eq. ( |T6"D is an analytic function in A in the vicinity of A = 0. It 
exhibits no cut or other singularities near the threshold and therefore does not contribute 
to the spectral density. We turn now to the second part II - (p). In contrast to the previous 
case, the integrand of this part contains H~(ipx) which behaves like a rising exponential 
function at large x, 

H-{ipx) ~ x ~ 1/2 e px . (18) 
Therefore the integral is represented by 

/•OO 

nx(Af-A)~/ x~ a e- Ax dx. (19) 

J A 

The function II~ (M — A) is non-analytic near A = because for A < the integrand 
in Eq. ([19]) grows in the large x region and the integral diverges at the upper limit. 
Therefore the function which is determined by this integral is singular at A < (E > 0) 
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and requires an interpretation for these values of the argument A. In the complex A 
plane with a cut along the negative axis the function is analytic. This cut corresponds 
to the physical positive energy cut. The discontinuity across the cut gives rise to the 
non- vanishing spectral density of the diagram. 

Let us first discuss the analytic part H + (p) of the diagram. This part reduces to a 
regular water melon diagram. Indeed, using the relation 



between Bessel functions of different kinds one can replace the Hankel function H^(ipx) 
with the McDonald function K\(px). Since the propagator of a massive particle (massive 
line in the diagram) is given by the McDonald function up to a power in x, this substitution 
shows that the weight function behaves like a propagator of an additional line with the 
"mass" p. The explicit expression is given by 



n + (a;) = U(x)D(x,p) is the polarization function of a new effective diagram which is equal 
to the initial polarization function multiplied by a propagator with p as mass parameter. 
We thus end up with a vacuum bubble of the water melon type with one additional line 
compared to the initial diagram (see Fig. |l|) . These diagrams have no singular behaviour 
at the production threshold p = M. As mentioned above, IT + (p) is analytic in A near 
the origin A = and can therefore be omitted in the calculation of the spectral density. 
All derivatives of II + (p) = II + (M — A) with respect to A are represented as vacuum 
bubbles with one additional line carrying rising indices. Such diagrams can be efficiently 
calculated within the recurrence relation technique developed in Ref. 

Therefore one is left with the second part IT~ (p) containing the non-analytical contri- 
butions in A which leads to a non-vanishing spectral density. Still the integral in Eq. (|T4]) 
cannot be done analytically. In order to obtain an expansion for the spectral density near 
the threshold in an analytical form we make use of the asymptotic series expansion for 



K x (z) = -e iX */ 2 H+(zz) 



(20) 




(21) 
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Figure 1: Representation of H + (p) as vacuum bubble with added line. The cross denotes 
an arbitrary number of derivatives to the specified line. 

the function U(x) which crucially simplifies the integrands but still preserves the singular 
structure of the integral in terms of the variable A. The asymptotic series expansion 
to order N of the main part of each propagator, i.e. of the McDonald function, is given 

by H 



7T \V2 

K as ( 



N ^(\,n) (A, TV) 
to(^) n (2*)" 



(A, n) := T ^ +n ~ 1/2 \ (22) 
v ' ; n\T(X-n- 1/2) v ; 



,1z. 

(9 G [0,1]). Therefore the asymptotic expansion of the function U(x) consists of an 
exponential factor e~ Mx and an inverse power series in x up to an order N which is 
closely related to N. It is this asymptotic expansion that determines the singularity 
structure of the integral. We write the whole integral in the form of the sum of two terms, 

n-(p) = v x+1 j (^pj X H^{ipx){Tl{x)-T^{x))x 2X+l+2£ dx (23) 

+tt a+1 J {^-\ A H^(ipx)Il™{x)x 2X+1+2£ dx = U di {p) + U as {p). 



The integrand of the first term (p) behaves as l/x N at large x while the integrand of 
the second term accumulates all lower powers of the large x expansion. Note that only the 
large x behaviour is essential for the near threshold expansion of the spectral density. This 



fact has been taken into account in Eqs. (pl|) and (|19"D where we introduced a cutoff A. 
However, from the practical point of view the calculation of the regularized integrals with 
an explicit cutoff is inconvenient. The final result of the calculation - the spectral density 
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of the diagram - is independent of the cutoff, but the integration is technically complicated 
if the cutoff is introduced. However, in extending the integration over the whole region of 
the variable x without using the cutoff one immediately encounters divergences at small 
x because the asymptotic expansion is invalid in the region near the origin, so one is 
not allowed to continue the integration to this region. The standard way to cope with 
such a situation is to introduce dimensional regularization. It allows one to deal with 
divergent expressions at intermediate stages of the calculation and is technically simple 
because it does not introduce any cutoff and therefore does not modify the integration 
region drastically. Note that dimensional regularization does not necessarily regularize all 
divergences in this case (in contrast to the standard case of ultraviolet divergences) but 
nevertheless suffices for our purposes. We use a parameter e to regularize the divergences 
at small x. This regularization prescription is an unorthodox version of the dimensional 
regularization (see e.g. Refs. |39[| ). 

The first part I\ dl {p) in Eq. ( ^3|) containing a difference of the polarization function 
and its asymptotic expansion as the integrand gives no contributions to the spectral den- 
sity up to a given order of the expansion in A. This is because the subtracted asymptotic 
series to order N cancels the inverse power behaviour of the integrand to this degree N. 
The integrand decreases sufficiently fast for large values of x and the integral converges 
even at A = 0. Therefore this term is inessential when the expansion of the spectral den- 
sity is evaluated up to some order. One can readily determine the order of the expansion 
near A = at which a contribution to the spectral density appears in using some further 
simplifications of the integrand of the term U dl (p) in Eq. (23). Namely, we replace the 
Hankel function under the integration sign by its asymptotic series expansion. The re- 
sulting exponential factor e^ _Af ^ can then be expanded in the parameter A = M — p and 
integrated together with the finite inverse power series in x. One obtains a finite power 
series in this parameter A which leads to a non-regular term of order A^ (for instance, 
A N In A or A N y/A). Therefore the part U di (p) is regular and gives no contribution to 
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the spectral density up to the order A . For this reason we concentrate on the expansion 
of the second part U as (p) and find that only this part contains the contribution to the 
spectral density up to the N. 

Therefore the expansion of the spectral density at small E is determined only by the 
integral Il as (p) of Eq. (^3|). This integral is still rather complicated to compute but we 
can go a step further in its analytical evaluation. Indeed, since the singular behaviour of 
U as (p) is determined by the behaviour at large x, we can replace the first factor, i.e. the 
Hankel function, in the large x region by its asymptotic expansion up to some order N. 
We use 



2 \ 1 / 2 

H7%(iz) = (— 1 e z+iXn/2 



(~i) jv (A, n) 

n=0 (2*)» (2s)" 



(24) 



(cf. Eq. (|22|) for the notation) to obtain a representation 



U das (p) = 7r A+1 |(^) H^(ipx)U a ^(x)x 2X+1+2e dx. (25) 

The index ll das n stands for "double asymptotic" and indicates that the integrand in 
Eq. (p5[) consists of a product of two asymptotic expansions: one for the polarization 
function H{x) and another for the Hankel function H^(x) as weight (or kernel). Both 
asymptotic expansions are straightforward and can be obtained from standard handbooks 
on Bessel functions. We therefore arrive at our final result: the integration necessary for 
evaluating the near threshold expansion of the water melon diagrams reduces to integrals 
of the type of Euler's Gamma function, i.e. integrals containing exponentials and powers. 
Indeed, the result of the expansion in Eq. ( p5|) is an exponential function e~ Ax times a 
power series in 1/x, namely 



JV-l 



a+2e„-Ax 



.4 



x e 



E -i (26) 



3 =o & 



where a has already been defined in Eq. ([H]) and the coefficients Aj are simple functions 
of the momentum p and the masses m r The expression in Eq. ( p6|) can be integrated 
analytically using 

fan 

x- a+2£ e- Ax dx = r(l - a + 2e)A a - 1 - 2e . (27) 
14 



The result is 

N-l 

U das (M - A) = AjT{l - a - j + 2e)A a+j - 1 - 2e . (28) 

j=0 

This expression is our final representation for the part of the polarization function of 
a water melon diagram necessary for the calculation of the spectral density near the 
production threshold. It is also one of the main results of our paper. 

Next we discuss the general structure of the expression in Eq. in detail. In the 
case where a takes integer values, these coefficients result in 1/e- divergences for small 
values of e. The powers of A in Eq. (|28|) have to be expanded to first order in e and give 

lA 2£ = l + lnA + 0(5). (29) 

Because of 

Disc ln(A) = ln(-£ - zO) - ln(-£ + iO) = -2m6{E) (30) 

II das (M — A) in Eq. ( p8|) contributes to the spectral density. For half-integer values of a 
the power of A itself has a cut even for e = 0. The discontinuity is then given by 



Disc = -2iVE0{E). (31) 

Our method to construct a threshold expansion thus simply reduces to the analytical 
calculation of the integral in Eq. (^) which can be done for arbitrary dimension and 
an arbitrary number of lines with different masses. In the next subsections we use our 
technique to work out some specific examples which demonstrate both the simplicity and 
efficiency of our method. 

3.1 Equal mass sunset diagram 

The polarization function represented by the sunset diagram with three propagators with 
equal masses m in D = 4 space-time dimensions is given by 

„. , m 3 Ki(mx) 3 , . 

n <*> = < 32 > 
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The exact spectral density is given by the integral representation in Eq. (§) which for this 
particular case reads 

2<7r rc+ioo 

P (s) = —= / / 1 (Cv / ^)n(C)C 2 rfC- (33) 

%\J S J c—ioo 

In order to obtain a threshold expansion of the spectral density in Eq. ( |33| ) we use Eq. (|28| ) 
to calculate the expansion of the appropriate part of the polarization function. To illus- 
trate the procedure we present the explicit shape of the integrand in Eq. ( p5|) which is 
given by an asymptotic expansion at large x, 

- 2 (f ) ^wdr = (4 ; )3p 3 /2 ^- 3+2£ >< 

x |l+ 9 3 + 9 27 15 +0(2T 3 )) (34) 

| 8ms 8px 128m 2 x 2 QAmpx 2 128p 2 x 2 j 

From Eq. (|34|) we can easily read off the coefficients Aj that enter the expansion in 
Eq. (p6|). The spectral density is obtained by performing the term-by-term integration of 
the series in Eq. ([34]) and by evaluating the discontinuity across the cut along the positive 
energy axis E > 0. The result reads 

£ 2 f n 1 7 2 3 3 39 4 57 5 
= ^V!! 1 - 2 V + 16 77 - 8^ + ^ - ^ (35) 
+ i^^ 6 _ ±j _ iOfL^s + 18603 _ 248829 1 
1024' 256' 32768' 65536' 524288' w ; J 

where the notation 77 = E/M, M = 3m is used. The simplicity of the derivation is 
striking. By no cost it can be generalized to any number of lines, arbitrary masses, and 
any space-time dimension. The standard equal mass sunset is chosen for the defmiteness 
only. It also allows us to compare our results with results available in the literature. 
Eq. (|35| ) reproduces the expansion coefficients dj obtained in Ref. |^TJ (the fourth column 
in Table 1 of Ref. [|3l[]) by a direct integration in momentum space within the technique 
of region separation |40| . 

The case of the equal mass standard sunset diagram is the simplest one. There exists 
an analytical expression for the spectral density of the sunset diagram with three equal 
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mass propagators in D = A space-time dimensions in terms of elliptic integrals [Ell] (see 



also Ref This expression can be used for a comparison with our exact result in 

Eq. (|33|) or with the expansion in Eq. ([35|). However, we only present the result for D = 2 
in order to keep the resulting expressions in a reasonably short form (cf. Ref. In 
D = 2 space-time dimensions the spectral density for a sunset diagram with equal masses 
m can be readily obtained. We just use the exact expression for the spectral density in the 
convolution representation @ and proceed towards n = 3 equal masses. The convolution 
function for two spectral densities in D = 2 dimensional space-time (A = 0) reads 

p{s; si; s 2 ) = , 1 (36) 

2nd (s — Sx — s 2 Y - AsiS 2 



The two spectral densities one has to convolute are the spectral density of a correlator 
with two equal masses and the spectral density of a single massive line. While the latter 
is given by p(s;m 2 ) = 8{s — m 2 ), the former can be obtained from Eq. (|36|) by inserting 



si = s 2 = m 2 , 



So the convolution leads to 

2. 2. 2^ 1 /"n '// 



p(s] m 2 ; m 2 ) = 1 (37) 

2-kJs(s - Am 2 ) 



p(s; m ;m ;m ) 



^ 2 J ^ 2 yj(s -m 2 - t) 2 - 4m 2 tyft(t - Am 2 ) 
1 AV-s-mf dt 



47r2 J *™ 2 ^/t(t - 4m 2 )((^i + m) 2 - t)((yfe - m) 2 - t) 
Now we use the relation (cf. Ref. f4l|j ) 



(38) 



K(k 2 ), (39) 



41 v/(*-to)(t-*i)(*2-t)(t3-t) y/(t3-ti)(t 2 -t ) 



k 2 = (fr-^-fr) (4Q) 

(t 3 -ti)fe-to) V ^ 



1 We thank A. Davydychev for bringing these papers to our attention. 
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with £ 3 > i 2 > t > t\ > to and the definition of the complete elliptic integral of the first 
kind 

w / 2 dtp ^ fir 1 2 



(remark the difference in the definition) for t = 0, t\ = 4m 2 , t 2 = (y/s — m) 2 , and 
^3 — (y/s + m ) 2 to perform the integration in Eq. (|38|). We obtain 

fc2 = ((vj ~ m ) 2 ~ 4m 2 )(y / i + m) 2 
" ((a/s + m) 2 - 4m 2 )( v / i - m) 2 

and finally end up with 

p(s; m 2 ; m 2 ; m 2 ) = K ^ \ (43) 

27r 2 (y / s — m)J (y/s + m) 2 — 4m 2 

Therefore the spectral density in terms of the energy E reads (see e.g. Ref. f4"4]| ) 

p(E) = 1 K(k 2 ), 

2vr 2 (2m + E)^(4m + E) 2 - 4m 2 

((4m + E) 2 -4m 2 )(2m + E) 2 ' v 7 

By expanding the elliptic integral in terms of the threshold parameter E one reproduces 
the threshold expansion in Eq. (|35|) . The result for D = 4 space-time dimensions is 
expressible by the elliptic integrals with some rational functions as factors that makes 
the result a bit longer. Note that the representation in Eq. flP|) is understood to be an 
analytical expression for the spectral density. However, it is a matter of taste whether 
the representation through the elliptic integrals as in Eq. ( |44"D is considered simpler (or in 
a more analytical form) than the integral representation in Eq. The only objection 
against the latter which one can find in the literature is that the Bessel functions are 



complicated (see e.g. Ref. [30]). But after more than a century of intensive investigation 
they are well-known and no more complicated than the square root of the fourth order 
polynomial which is used in Eq. ( |39|) to define the elliptic integral. 
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3.2 Equal mass water melon diagrams 
with four or more propagators 

The water melon diagrams with four or more propagators cannot be easily done by using 
the momentum space technique because it requires the multiloop integration of entangled 
momenta. Within the configuration space technique the generalization to any number of 
lines (or loops) is immediate by no effort. Consider first a three- loop case of water melon 
diagrams (also called banana diagrams |28] or basketball diagrams ]nj). The polarization 



function of the equal mass water melon diagram with four propagators in D = 4 space- 
time is given by 

n W = !=W ,45) 



The exact spectral density of this diagram can be obtained from Eq. (|33| ) while the near 
threshold expansion can be found using Eq. (|28|) . We construct the expansion of the 
spectral density near threshold explicitly and compare it with the exact result. The 
expansion of the integrand (cf. Eq. (p5j)) reads 

f " 3 3 3 15 9 0(H- (46) 



[ 2mx 8px 8m 2 x 2 I28p 2 x 2 lQmpx 2 
After the integration and the calculation of the discontinuity we obtain the expansion of 
the spectral density in the form 

E 7 / 2 M 1 ' 2 r 1 81 2 2811 3 17581 4 

p(E) = = 1 1 77 H n 2 7T 3 H 77 4 (47) 

268807T 5 v / 2l 4 ' 352 ' 18304 ' 292864 ' v ' 

1085791 . 597243189 fi 4581732455 7 496039631453 R nl q ,l 
-rf rf H rf r/ 8 + 0(rf) ) 



19914752 ' 3027042304 ' 12108169216 ' 810146594816 
where r] — E/M and M = 4m is the threshold value. One sees the difference with the 
previous three-line case. In Eq. (f47|) the cut represents the square root branch while in 
the three-line case it was a logarithmic cut. One can easily figure out the reason for this 
by looking at the asymptotic structure of the integrand. For even number of lines (i.e. 
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odd number of loops) it is a square root branch, while for an odd number of lines (even 
number of loops) it is a logarithmic branch. This is true in even space-time dimensions. In 
the general case the structure of the cut depends on the dimensionality of the space-time 
as well. The general formula reads 

p{E) ~ £(A+i/2)(«-i)-i(i + 0(E)). (48) 

For D = 4 space-time dimension (i.e. A = 1) we can verify the result of Ref. ]7| (cf. 
Eq. @D), 

p(E) ~£ (3n - 5)/2 (l+0(£)). (49) 
Numerically Eq. fl4"?D reads 

p(E) = 8.5962- 10" 5 £ 7/2 M 1/2 {l.000- 0.250// + 0.230r? 2 (50) 
— 0.154?7 3 + 0.060r/ 4 + 0.055r? 5 - 0.197^ + 0.378r? 7 - 0.612^ + 0(r/ 9 )} 

where we have written down the coefficients up to three decimal places. It is difficult 
to say anything definite about the convergence of this series. By construction it is an 
asymptotic series. However, we stress that the practical (or explicit) convergence can 
always be checked by comparing series expansions like the one shown in Eq. ([5(]) with the 
exact spectral density given in Eq. (|33| ) by numerical integration. 

We conclude this part of the paper by the statement that the spectral density of the wa- 
ter melon diagram is most efficiently calculable within the configuration space technique. 
Whether it is the exact result or the expansion, the configuration space technique can 
readily deliver the desired result. The exact formula in Eq. (^) as well as the threshold 
expansion obtained from it can be used to calculate the spectral density for an arbitrarily 
large number of internal lines. We stress that the case of different masses does not lead to 
any complications within the configuration space technique: the exact formula in Eq. (|8]) 
and/or the near threshold expansion work equally well for any arrangement of masses. 
We do not present plots for general cases of different masses because they are not very 
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illustrative, showing only the common threshold. However, there is some interesting kine- 
matic regime for different masses which is important for applications and which, to our 
best knowledge, have not been touched earlier. An analytical solution for the expansion 
of the spectral density in this regime is given in the next section. 




0.2 0.4 0.6 0.8 1 

7] = E/M 

Figure 2: Various results for the spectral density for n = 3 equal masses in D = 4 space- 
time dimensions in dependence on the threshold parameter E/M. Shown are the exact 
solution obtained by using Eq. ( |3"3"D (solid curve) and threshold expansions for different 
orders taken from Eq. (^) (dashed to dotted curves). 



4 Strongly asymmetric case mo <C M 

The threshold expansion for equal (or close) masses breaks down for E « M = J2 m i- 
The example is shown in Fig. |2| for the D = 4 proper sunset. However, if the masses are 
not equal, the region of the break-down of the expansion is determined by the mass with 
the smallest numerical value. The simplest example where one can see this phenomenon 
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is the analytical expression for the spectral density of the simple loop (degenerate water 
melon diagram) with two different masses m\ and m 2 . In D = 4 space-time dimensions 
(see e.g. Ref. J7|) one has 



. /E(E + 2m 1 )(E + 2m 2 )(E + 2M) 

P{E) = (An(M + E)) 2 (51) 

where M = mi + m-i- The threshold expansion is obtained by expanding the right hand 
side of Eq. (|5T]) in E for small values of E. If m 2 is much smaller than mi, the expansion 
breaks down at E ~ 2m 2 . The break-down of the series expansion can also be observed 
in more general cases. If one of the masses (which we call mo) is much smaller than the 
other masses, the threshold expansion is only valid in a very limited region E ^ 2mo- 

To generalize the expansion and extend it to the region of E ~ M one has to treat the 
smallest mass exactly. In this case one can use a method which we call the resummation 
of the smallest mass contributions. Below we describe the resummation technique. We 
start with the representation 

IF™(p) = n x+1 J X H^(ipx)UZ (x)x 2X+1 ^dx (52) 

which is the part of the polarization function contributing to the spectral density. The 
integrand in Eq. ( |52|) has the form 

U% {x)=U%_ 1 (x)D{m Qt x) (53) 

where the asymptotic expansions are substituted for all the propagators except for the 
one with the small mass rriQ. This is indicated by the index u pas" in Eq. ( |52"D which stands 
for "partial asymptotic". The main technical observation leading to the generalization 
of the expansion method is that U pas (p) is still analytically computable in a closed form. 
Indeed, the genuine integral to compute has the form 

oo 

x»- l e- &x K u {f3x)dx = 

(2a)^ I> + 1/2) 21 l^"' 2 ' " + 2 ' 1 ~ (M) 
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where a = A — mo and (3 = mo. The integral U pas (p) in Eq. (|52| ) is thus expressible in 



terms of hypergeometric functions [ 25 , W\. For constructing the spectral density, being our 
main concern as mentioned before, one has to find the discontinuity of the right hand side 
of Eq. flSip . There are several ways to do this. We proceed by applying the discontinuity 
operation to the integrand of the integral representation of the hypergeometric function. 
The resulting integrals are calculated again in terms of hypergeometric functions. Indeed, 



—Disc / x^e^KJfix^x 
2m Jo 



r(3/2) ( l-y-v 2-n-v 3 /3 2 \ 

«V2-^ r(3/2 - n) 2 ^ { 2 ' 2 '2 ^a*) ^ 

where a = E + m . The final expression in Eq. (FJ^) completely solves the problem of the 
generalization of the near threshold expansion technique. For integer values of fi there 
are no singular Gamma functions (with negative integer argument). Therefore we can lift 
up the regularization and set e = when using this expression. We thus have found a 
direct transition from the polarization function as expressed through the integral to the 
spectral density in terms of one hypergeometric function for each genuine integral. There 
is no need to use the recurrence relations available for hypergeometric functions. 

In the following subsections we give explicit examples for D = 4 and compare with the 
exact result in Eq. (|33D and the pure expansion near the threshold. In the following the 
standard threshold expansion without resummation is called the pure threshold expansion. 



4.1 The two-line water melon with a small mass 

We start with a (over) simplified example of the two- line diagram with masses m and 
rrio C m in four space-time dimensions. In this example the expansion of the spectral 
density and its generalized expansion can be readily compared analytically with the exact 
result in Eq. (|5T|). This is the feature that justifies our discussion in this section. The 
results for the spectral density of this diagram are shown in Fig. |3|. The solid curve displays 
the exact result obtained by using Eq. ([33]) (which reproduces the analytical expression 
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Figure 3: Various solutions for the spectral density for two masses m and mo <C m and 
= 4 space-time dimensions. Shown are the exact solution which is obtained by using 
Eq. ( |3~3"1) (solid curve), the pure threshold expansions using Eq. (|^) (dotted curves), and 



the solutions for the resummation of the smallest mass contributions like in Eq. (|57|) 
(dashed curves), both expansions from the first up to the fourth order in the asymptotic 
expansion. For the pure threshold expansion the order is indicated explicitly. 

in Eq. (|51|)). We compare this result with the two expansions. 

The pure expansion of the spectral density near threshold (the second order asymptotic 
expansion should suffice to show the general features in a short and concise form) is given 
by 



p das ( 







{• 








1 


12 




+ 


m 2 


m m 



m 



E 
1 



(56) 



where M = m + m . As mentioned above, this series breaks down for E > 2m (see 
Eq. (HD). If we look at the dotted curves in Fig. |] this becomes obvious. Here we have 
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plotted the series expansions up to the fourth order with the mass arrangement mo = 
m/10. The dashed lines represent the resummation of the smallest mass contributions. 
The analytical expression for the spectral density of the polarization function in Eq. (|52|) 
for the generalized asymptotic expansion based on Eq. ([55]) is given by 



>mE(E + 2m ) i / 13 m 2 \ 

• {E) ~ 8n 2 (E + Mf/ 2 ( 2 1 \ , 2'2' (E + m ) 2 ) 

E(E + 2mo) (1 5 m 2 \ 

+ 8m(£ + M) 2 1 \j2' ' 2' (E + mo) 2 ) 1 J 

£ 2 (£ + 2m ) 2 / , 16m(E + M)\ ( 3 7 1 \ l 

1 + r/T-i , 2-ri i, -; -; J- 



128m 2 (E + M) 2 V 5(£ + m ) 2 / V 2 2 (E + m ) 2 J ' J 

We have set the regularization parameter e = because the spectral density is finite. With 
e = the resulting expressions for the hypergeometric functions in Eq. Q55D simplify. The 
first term in the curly brackets of Eq. is obviously equal to 1 in this limit because 
the first parameter of the hypergeometric function vanishes for e — 0. However, we keep 
Eq. (|57D in its given form to show the structure of the contributions. The generalized 
threshold expansion has the form 

f? as (E) = g (E, mo) + E 9l (E, m ) + E 2 g 2 (E, m ) + . . . (58) 

where the functions gj{E,mo) represent effects of the resummation of the smallest mass 
and are not polynomials in the threshold parameter parameter E. In the simple two-line 
case the hypergeometric functions reduce to elementary functions. For instance, 

/ 1 5 m„ \ r \ 

= (59) 

3(E + m ) ( m 2 , / ^ + m + J E(E + 2m ) ' 

1 E + m In 



2£(£ + 2m ) ^ 2^/E(E + 2m ) \^ + m o - ^E(E + 2m ) 

Higher order contributions are given by hypergeometric functions with larger numerical 
values of the parameters. They can be simplified by using Gaussian recurrence relations 
for hypergeometric functions (see e.g. Ref. 
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resummed expansions 
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0.5 1 1.5 2 

rj = E/ (m + m ) 

Figure 4: The same as Fig. |3| where the spectral density is normalized to the leading order 
expression of the pure threshold expansion. 

The convergence of the expansion in Eq. ([57]) breaks down only at E ~ M = m + itlq. 
The resummation leads to an essential improvement of the convergence in comparison 
with the pure threshold expansion. In Fig. [| we show the same curves divided by the 
leading order term. This representation is more convenient for the diagrams which we 
will discuss in following subsections. 

Note that Eq. (|5^) does not lead to the exact function in Eq. ([51]) because terms of 
order E N stemming from the difference part I\ dl {p) of the correlator are missing. It simply 
corrects the behaviour of the coefficient functions by the small mass contributions. 

4.2 The sunset diagram with a small mass 

Here we analyze the sunset diagram with two equal masses m and a third mass mo <C m 
(m = m/10). The exact result obtained by using Eq. ( |55p and normalized to the leading 
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Figure 5: The spectral density for the sunset diagram in D = 4 space-time dimensions with 
a tiny mass mo, normalized to the general power behaviour. Shown are the exact result 
obtained by using Eq. (|3~3"l) (solid curve), the threshold expansion according to Eq. (|60|) 
(dotted curves), and the result for the resummation of the smallest mass contributions 
according to Eq. (|BTD (dashed curves). 

order term is shown in Fig. [5] as the solid curve. The pure expansion near the threshold 
reads 



~ das(F] mEW^M( /J_ I _ 13\E 
P 1 J 128tt 3 M 2 \ U m M) 8 



m m 
5 4 39 153 



1115\ E 2 



mo 



+ 0(E 3 



(60) 



l mom ' mM M 2 J 512 j 

It is shown by the dotted curves in Fig. [5[ In case of the resummation of the smallest 
mass contributions we obtain hypergeometric functions which do not obviously reduce to 
elementary functions in this case. The result for the spectral density within the asymptotic 
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expansion up to the second order in Eq. ( p2] ) is given by 

mE 2 (E + 2m ) 2 r 



[P as (E) 



mS 



5i2tt 3 (£ + m ) 3 / 2 (£ + M) 3 / 2 | x \ 4' 4 



+ 



3m 



g(g + 2m ) / 

8m(E + M) I ' 2(£ + m ) 



4' 4' 



(£ + m ) 2 



1 + 



5m 



£ 2 (£ + 2m ) 2 

512m 2 (£ + M) 2 y ' 2(E + m 



9m 



2(E + m ) 



x 



^(j,j;5;i- 



(61) 



(F + m ) 2 ; J 

We see that the dashed curves in Fig. ^ that represent the result in Eq. ( pip approximate 
the exact curve much better than the dotted curves. 
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Figure 6: The spectral density for the four-line water melon diagram in D = 4 space-time 
dimensions with a tiny mass mo, normalized to the general power behaviour. Shown 
are the exact result obtained by using Eq. ( |3~3D (solid curve), the threshold expansion 
according to Eq. (|62|) (dotted curves), and the result for the resummation of the smallest 
mass contributions according to Eq. ( |63D (dashed curves). 
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4.3 The four-line water melon with a small mass 



With this example we conclude our consideration of the strongly asymmetric case and at 
the same time show the way to the multi-line water melon diagrams which can be treated 
in an analogous manner. The result for the exact expression obtained by using Eq. (|33|) 
is shown in Fig. ^ as a solid line, normalized to the leading order term. The dotted lines 
represent the results for the pure expansion near threshold which is given by 



m 



3/ 



/ 2 £ 7 / 2 V2mo" L ( 1 3 19 \ £ 



5 3 28 368 2195\ E 2 



\mQ m 2 mom mM M 2 J 1056 
The asymptotic expansion to the second order in Eq. fl5"2|) gives 



0(E 3 



26880vr 5 (£ + m ) 3 (£ + M) 3 / 2 j 2 ^2' '2' (E + m ) 2 
E(E + 2m ) f 1 + 8m \ t 5 11 m 2 \ 



8m(E + M)\ 3(E + m )J \ ' 2' 2 ' (E + m ) 2 ) 
E 2 (E + 2m y ( 32m 2 \ (5 13 m 2 \\ 

^U08(E + M) 2 \ 3(E + m ) 2 ) 2 1 ^2' ' 2 ' (E + m ) 2 ) J" 1 J 

In Fig. |^ one can see how the expansion improves if the resummation of the smallest mass 
contributions (displayed as dashed lines) is performed. 

The result of this section is quite general and applicable to all cases of one small mass. 
For some particular arrangement of masses one can obtain even simpler expressions as 
discussed in the next section. 



4.4 The convolution with a small mass 

In this section we obtain a result for the resummation of the smallest mass effects along 
a different route, namely, via the convolution of spectral densities. However, this method 
works in a narrower kinematic region than the method described in the previous section. 
In D = 4 space-time dimensions, the convolution weight is given by 

p(s; Si, s 2 ) = (4^2^ \A S ~ Sl ~ S2 > 2 ~ 4:SlS2 - ( 64 ) 
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The upper limit of the integration is determined by the requirement of positivity of the 
the square root argument. The zeros of the square root with respect to s 2 ar e given by 
s f = (y/s =t V^i) 2 > an d the demand (s 2 — s^)(s2 — s 2 ) > together with s 2 > s 2 leads 
to s 2 > S2 or s 2 < s 2 • The physical region is the latter one. With pi(s) = S(s — TOq) for 
the spectral density of the single small mass line we obtain 



where the low limit of integration is M' = M — m . We insert s = (M + E) 2 and 
s 2 = (M' + E'f and obtain 



where p'(E') = p 2 ((M' + E') 2 ). For this function we use the threshold expansion in E' /M' 
as expansion parameter. For small E < M' the threshold expansion inserted for p'(E') is 
valid because E < M' implies E' < M'. The described procedure can be extended to the 
case of a very light sub-block of the diagram, e.g. a light fish diagram. In this case we 
have to replace pi(s) by the spectral density of the light sub-diagram which is well-known. 



The analytic structure of water melon diagrams is completely fixed by the dispersion 
representation. Therefore we have focussed on the computation of the spectral density as 
the basic quantity important both for applications and the theoretical investigation of the 
diagram. However, with an analytical expression for the spectral density p(s) at hand we 
can readily reconstruct the non-analytic piece of the polarization function in momentum 




(65) 




(66) 



5 Recovering U(p) through p(s) near threshold 
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space by using the dispersion relation 

nM = /^. (67) 

J s — p z 

We rewrite this equation in terms of threshold parameters according to p = M — A and 
s = (M + E) 2 and obtain 

^t/ax A x Z" 00 2(M + E)p(E)dE . , 

II A = n M - A = / i — x , 7 -. 68 

v ; v ; io (.B + A)(2M + -E — A) v ; 

UV singularities can be removed by subtraction or by dimensional regularization. We 
again use the unorthodox dimensional regularization prescription. For a general form of 
the threshold expansion p(E) = E 1 a>kE k we have to calculate integrals of the form 

fr ( M = f 00 E ° dE vr A- - (2M - A)- 

1 J " Jo (E + A)(2M + E — A) sin(Trcr) 2(M — A) { } 

Only the powers A CT contribute to the singular part of the polarization function. Ex- 
pressions like the one presented in Eq. (|69|) then allow one to restore that part of the 
polarization function H(p) which has singularities near the threshold. 



6 Conclusion 

We have discussed the configuration space technique for the calculation of ra-line two- 
point diagrams termed water melon diagrams. This technique allows one to calculate the 
spectral density for arbitrary space-time dimensions and any number of internal lines with 
arbitrary mass values. Within this technique one can use either an exact representation 
as one-fold integral or an expansion near the production threshold. We have developed an 
efficient method for constructing the near threshold expansion of water melon diagrams 
that uses only asymptotic expansions of Bessel functions which are well-known and simple 
functions. We have considered a strongly asymmetric mass arrangement where one mass 
is much smaller than the others. In this case the "pure" threshold expansion which is 
done in terms of the threshold parameter E breaks down at the energies in the vicinity 
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Figure 7: The spectral density for the four-line water melon diagram with equal masses 
for D = e = 2.718 . . ., D — 3, and D = n = 3.14 . . . space-time dimensions. 

of the smallest mass. In order to extend the approximation to higher energies we have 
introduced a generalized threshold expansion which exactly resums all contributions of 
the smallest mass. We have presented closed expressions for the generalized expansion in 
several particular cases and demonstrated the improvement of the convergence gained by 
the resummation of the smallest mass contributions. The particular kinematic regime of 
this case could be treated analytically because it reduced to the evaluation of the one-fold 
integral in terms of hypergeometric functions. The discontinuity across the physical cut 
for the generalized expansion has been found in terms of hypergeometric functions as well. 

To conclude, we stress that the configuration space technique is a powerful and con- 
venient tool for investigating different properties of water melon diagrams. The practical 
convenience of our method is demonstrated in Fig. [7| where we have plotted the spec- 
tral density for a four-line water melon diagram in D = e = 2.718..., D — 3, and 
D — 7T = 3.14 . . . space-time dimensions. 



32 



Acknowledgments 

We thank Andrei Davydychev and Jiirgen Korner for interesting discussions and care- 
ful reading the manuscript. A. A. Pivovarov acknowledges a valuable discussion with 
V.A. Smirnov on the region separation technique in threshold expansions. The work 
of A. A. Pivovarov is supported in part by the Volkswagen Foundation under contract 
No. 1/73611 and by the Russian Fund for Basic Research under contract 99-01-00091. 
S. Groote gratefully acknowledges a grant given by the Graduiertenkolleg "Eichtheorien 
- experimented Tests und theoretische Grundlagen" . 

References 

[1] Particle Data Group, "Review of Particle Properties", Eur. Phys. J. C3 (1998) 1 

[2] K.G. Chetyrkin, J.H. Kiihn and A. Kwiatkowski, Phys. Rep. 277 (1996) 189 

[3] K.G. Chetyrkin, J.H. Kiihn and M. Steinhauser, Nucl. Phys. B505 (1997) 213 

[4] D.J. Broadhurst, Eur. Phys. J. C8 (1999) 311 

[5] O.M. Ogreid, P. Osland, "More series related to the Euler series", 
Report No. BERGEN-1999-04, |hep-th/ 9904206 



[6] R. Coquereaux, Phys. Rev. D23 (1981) 2276; 
V.I. Zakharov, Nucl. Phys. B385 (1992) 452; 
D.J. Broadhurst and A.G. Grozin, Phys. Rev. D52 (1995) 4082 

[7] S. Groote, J.G. Kdrner and A.A. Pivovarov, Nucl. Phys. B542 (1999) 515; 
Eur. Phys. J. Cll (1999) 279; Phys. Lett. 443 B (1998) 269 

[8] F.A. Berends, M. Buza, M. Bdhm and R. Scharf, Z. Phys. C63 (1994) 227 

[9] F.A. Berends, A.I. Davydychev and N.I. Ussyukina, Phys. Lett. 426 B (1998) 95 

33 



[10] S. Coleman and E. Weinberg, Phys. Rev. D7 (1973) 1888; 
R. Jackiw, Phys. Rev. D9 (1974) 1686; 
R. Jackiw and S.Templeton, Phys. Rev. D23 (1981) 2291 

[11] A.K. Rajantie, Nucl. Phys. B480 (1996) 729 

[12] J.M. Chung and B.K. Chung, "Renormalization group improvement of the effective 
potential in massive $ 4 theory: NNNLO logarithm resummation" , Report No. MIT- 
CTP-2929, |hep-th/991li96| 



[13 

[14 
[15 

[16 
[17 

[18 



J.O. Andersen, E. Braaten and M. Strickland, "The massive thermal basketball dia- 
gram", [Me^pTi/0002048| 



J.M. Chung and B.K. Chung, Phys. Rev. D59 (1999) 105014 

D.J. Gross, R.D. Pisarski and L.G. Yaffe, Rev. Mod. Phys. 53 (1981) 43; 
T. Appelquist and R.D. Pisarski, Phys. Rev. D23 (1981) 2305 

B.L. Ioffe, Nucl. Phys. B188 (1981) 317 

N.V. Krasnikov, A. A. Pivovarov and A.N. Tavkhelidze, 
Z. Phys. C19 (1983) 301; JETP Lett. 36 (1982) 333 

A. A. Ovchinnikov, A. A. Pivovarov and L.R. Surguladze, 
Yad. Fiz. 48 (1988) 562; Int. J. Mod. Phys. A6 (1991) 2025 



[19] S. Groote, J. G. Korner and A. A. Pivovarov, Phys. Rev. D61 (2000) 071501 

[20] G. 't Hooft, Nucl. Phys. B72 (1974) 461 

[21] E. Witten, Nucl. Phys. B160 (1979) 57 

[22] S.A. Larin et al, Yad. Fiz. 44 (1986) 1066 

[23] S.A. Kulagin and A.A. Pivovarov, Yad. Fiz. 45 (1987) 952 

34 



[24 
[25 
[26 
[27; 
[28 
[29 



[30 

[31 
[32 

[33 
[34 
[35 
[36 



P. Post and K. Schilcher, Phys. Rev. Lett. 79 (1997) 4088 

J. Gasser and M.E. Sainio, Eur. Phys. J. C6 (1999) 297 

S. Narison and A.A. Pivovarov, Phys. Lett. 327 B (1994) 341 

I.I. Balitsky, D.I. D'Yakonov and A.V. Yung, Phys. Lett. 112 B (1982) 71 

D.J. Broadhurst, Z. Phys. C54 (1992) 599 

P. A. Baikov, "Advanced method of solving recurrence relations for multiloop Feyn- 
man integrals" , talk presented at the 6 th International Workshop on New Computing 
Techniques in Physics Research (AIHENP 99), Heraklion, Greece, 12-16 April, 1999 

N.E. Ligterink, "Solving multiloop Feynman diagrams using light front coordinates", 
Report No. ECT-99-016, |hep-ph/99114TTl 

A. I. Davydychev and VA. Smirnov, Nucl. Phys. B554 (1999) 391 

B. Kastening, H. Kleinert, "Efficient algorithm for perturbative calculation of multi- 
loop Feynman integrals", |quant-ph/ 99090 17 



G.N. Watson, "Theory of Bessel functions", Cambridge, 1944 

E. Mendels, Nuovo Cim. 45 A (1978) 87 

C.S. Meijer, Proc. Amsterdam Akad. Wet. (1940) 599; 702 

A. Erdelyi (Ed.), "Tables of integral transformations", 
Volume 2, Bateman manuscript project, 1954 



[37] A.H. Hoang et al, "Top-antitop pair production close to threshold: Synopsis of recent 
NNLO results," Report No. SLAC-PUB-8369, frep-ph/0001286 ; 



A. Czarnecky and K. Melnikov, Phys. Rev. Lett. 80 (1998) 2531; 

M. Beneke, A. Signer and V.A. Smirnov, Phys. Rev. Lett. 80 (1998) 2535; 

35 



K. Melnikov and A. Yelkhovsky, Nucl. Phys. B528 (1998) 59; 

A.A. Penin and A. A. Pivovarov, |hep-ph/ 99042781 ; Nucl. Phys. B549 (1999) 217; 

Nucl. Phys. B550 (1999) 375; Phys. Lett. 443 B (1998) 264; 

O.I. Yakovlev, Phys. Lett. 457 B (1999) 170; 

A.H. Hoang and T. Teubner, Phys. Rev. D60 (1999) 114027; 

T. Nagano, A. Ota and Y. Sumino, Phys. Rev. D60 (1999) 114014 

[38] N.N. Bogoliubov and D. V. Shirkov, "Quantum fields" , Benjamin, 1983 

[39] A.A. Pivovarov, Phys. Lett. 236 B (1990) 214; Phys. Lett. 263 B (1991) 282 

[40] V.A. Smirnov, Phys. Lett. B404 (1997) 101; 

M. Beneke and V.A. Smirnov, Nucl. Phys. B522 (1998) 321; 
V.A. Smirnov and E.R. Rakhmetov, Teor. Mat. Fiz. 120 (1999) 64 

[41] A. P. Prudnikov, Yu.A. Brychkov and O.I. Marichev, 

"Integrals and Series", Vol. 2, Gordon and Breach, New York, 1990 

[42] B.Almgren, Arkiv for Fysik, 38 (1967) 161; 

S.Bauberger, F.A. Berends, M. Bohm and M. Buza, Nucl. Phys. B434 (1995) 383 

[43] S. Groote, J. G. Korner and A. A. Pivovarov, Phys. Rev. D60 (1999) 061701 

[44] A.A. Pivovarov and V.F. Tokarev, Yad. Fiz. 41 (1985) 524 

[45] I.S. Gradshteyn and I.M. Ryzhik, 

"Tables of integrals, series, and products", Academic Press, 1994 

[46] M. Abramowitz, LA. Stegun (eds.), "Handbook of Mathematical Functions", Dover 
Publ. Inc., New York, 9th Printing, 1970 



36 



